Proteogenomic characterization of difficult-to-treat breast cancer with tumor cells enriched through laser microdissection.

BACKGROUND: Breast cancer (BC) is the most commonly diagnosed cancer and the leading cause of cancer death among women globally. Despite advances, there is considerable variation in clinical outcomes for patients with non-luminal A tumors, classified as difficult-to-treat breast cancers (DTBC). This study aims to delineate the proteogenomic landscape of DTBC tumors compared to luminal A (LumA) tumors. METHODS: We retrospectively collected a total of 117 untreated primary breast tumor specimens, focusing on DTBC subtypes. Breast tumors were processed by laser microdissection (LMD) to enrich tumor cells. DNA, RNA, and protein were simultaneously extracted from each tumor preparation, followed by whole genome sequencing, paired-end RNA sequencing, global proteomics and phosphoproteomics. Differential feature analysis, pathway analysis and survival analysis were performed to better understand DTBC and investigate biomarkers. RESULTS: We observed distinct variations in gene mutations, structural variations, and chromosomal alterations between DTBC and LumA breast tumors. DTBC tumors predominantly had more mutations in TP53, PLXNB3, Zinc finger genes, and fewer mutations in SDC2, CDH1, PIK3CA, SVIL, and PTEN. Notably, Cytoband 1q21, which contains numerous cell proliferation-related genes, was significantly amplified in the DTBC tumors. LMD successfully minimized stromal components and increased RNA-protein concordance, as evidenced by stromal score comparisons and proteomic analysis. Distinct DTBC and LumA-enriched clusters were observed by proteomic and phosphoproteomic clustering analysis, some with survival differences. Phosphoproteomics identified two distinct phosphoproteomic profiles for high relapse-risk and low relapse-risk basal-like tumors, involving several genes known to be associated with breast cancer oncogenesis and progression, including KIAA1522, DCK, FOXO3, MYO9B, ARID1A, EPRS, ZC3HAV1, and RBM14. Lastly, an integrated pathway analysis of multi-omics data highlighted a robust enrichment of proliferation pathways in DTBC tumors. CONCLUSIONS: This study provides an integrated proteogenomic characterization of DTBC vs LumA with tumor cells enriched through laser microdissection. We identified many common features of DTBC tumors and the phosphopeptides that could serve as potential biomarkers for high/low relapse-risk basal-like BC and possibly guide treatment selections.

Mapping microhabitats of lignocellulose decomposition by a microbial consortium.

The leaf-cutter ant fungal garden ecosystem is a naturally evolved model system for efficient plant biomass degradation. Degradation processes mediated by the symbiotic fungus Leucoagaricus gongylophorus are difficult to characterize due to dynamic metabolisms and spatial complexity of the system. Herein, we performed microscale imaging across 12-µm-thick adjacent sections of Atta cephalotes fungal gardens and applied a metabolome-informed proteome imaging approach to map lignin degradation. This approach combines two spatial multiomics mass spectrometry modalities that enabled us to visualize colocalized metabolites and proteins across and through the fungal garden. Spatially profiled metabolites revealed an accumulation of lignin-related products, outlining morphologically unique lignin microhabitats. Metaproteomic analyses of these microhabitats revealed carbohydrate-degrading enzymes, indicating a prominent fungal role in lignocellulose decomposition. Integration of metabolome-informed proteome imaging data provides a comprehensive view of underlying biological pathways to inform our understanding of metabolic fungal pathways in plant matter degradation within the micrometer-scale environment.

Deciphering ApoE Genotype-Driven Proteomic and Lipidomic Alterations in Alzheimer's Disease Across Distinct Brain Regions.

Alzheimer's disease (AD) is a neurodegenerative disease with a complex etiology influenced by confounding factors such as genetic polymorphisms, age, sex, and race. Traditionally, AD research has not prioritized these influences, resulting in dramatically skewed cohorts such as three times the number of Apolipoprotein E (APOE) ε4-allele carriers in AD relative to healthy cohorts. Thus, the resulting molecular changes in AD have previously been complicated by the influence of apolipoprotein E disparities. To explore how apolipoprotein E polymorphism influences AD progression, 62 post-mortem patients consisting of 33 AD and 29 controls (Ctrl) were studied to balance the number of ε4-allele carriers and facilitate a molecular comparison of the apolipoprotein E genotype. Lipid and protein perturbations were assessed across AD diagnosed brains compared to Ctrl brains, ε4 allele carriers (APOE4+ for those carrying 1 or 2 ε4s and APOE4- for non-ε4 carriers), and differences in ε3ε3 and ε3ε4 Ctrl brains across two brain regions (frontal cortex (FCX) and cerebellum (CBM)). The region-specific influences of apolipoprotein E on AD mechanisms showcased mitochondrial dysfunction and cell proteostasis at the core of AD pathophysiology in the post-mortem brains, indicating these two processes may be influenced by genotypic differences and brain morphology.

Long-read RNA-seq atlas of novel microglia isoforms elucidates disease-associated genetic regulation of splicing.

Microglia, the innate immune cells of the central nervous system, have been genetically implicated in multiple neurodegenerative diseases. We previously mapped the genetic regulation of gene expression and mRNA splicing in human microglia, identifying several loci where common genetic variants in microglia-specific regulatory elements explain disease risk loci identified by GWAS. However, identifying genetic effects on splicing has been challenging due to the use of short sequencing reads to identify causal isoforms. Here we present the isoform-centric microglia genomic atlas (isoMiGA) which leverages the power of long-read RNA-seq to identify 35,879 novel microglia isoforms. We show that the novel microglia isoforms are involved in stimulation response and brain region specificity. We then quantified the expression of both known and novel isoforms in a multi-ethnic meta-analysis of 555 human microglia short-read RNA-seq samples from 391 donors, the largest to date, and found associations with genetic risk loci in Alzheimer's disease and Parkinson's disease. We nominate several loci that may act through complex changes in isoform and splice site usage.

Human cerebrospinal fluid contains diverse lipoprotein subspecies enriched in proteins implicated in central nervous system health.

Lipoproteins in cerebrospinal fluid (CSF) of the central nervous system (CNS) resemble plasma high-density lipoproteins (HDLs), which are a compositionally and structurally diverse spectrum of nanoparticles with pleiotropic functionality. Whether CSF lipoproteins (CSF-Lps) exhibit similar heterogeneity is poorly understood because they are present at 100-fold lower concentrations than plasma HDL. To investigate the diversity of CSF-Lps, we developed a sensitive fluorescent technology to characterize lipoprotein subspecies in small volumes of human CSF. We identified 10 distinctly sized populations of CSF-Lps, most of which were larger than plasma HDL. Mass spectrometric analysis identified 303 proteins across the populations, over half of which have not been reported in plasma HDL. Computational analysis revealed that CSF-Lps are enriched in proteins important for wound healing, inflammation, immune response, and both neuron generation and development. Network analysis indicated that different subpopulations of CSF-Lps contain unique combinations of these proteins. Our study demonstrates that CSF-Lp subspecies likely exist that contain compositional signatures related to CNS health.

TopPICR: A Companion R Package for Top-Down Proteomics Data Analysis.

Top-down proteomics is the analysis of proteins in their intact form without proteolysis, thus preserving valuable information about post-translational modifications, isoforms, and proteolytic processing. However, it is still a developing field due to limitations in the instrumentation, difficulties with the interpretation of complex mass spectra, and a lack of well-established quantification approaches. TopPIC is one of the popular tools for proteoform identification. We extended its capabilities into label-free proteoform quantification by developing a companion R package (TopPICR). Key steps in the TopPICR pipeline include filtering identifications, inferring a minimal set of protein accessions explaining the observed sequences, aligning retention times, recalibrating measured masses, clustering features across data sets, and finally compiling feature intensities using the match-between-runs approach. The output of the pipeline is an MSnSet object which makes downstream data analysis seamlessly compatible with packages from the Bioconductor project. It also provides the capability for visualizing proteoforms within the context of the parent protein sequence. The functionality of TopPICR is demonstrated on top-down LC-MS/MS data sets of 10 human-in-mouse xenografts of luminal and basal breast tumor samples.

Mass spectrometry-based direct detection of multiple types of protein thiol modifications in pancreatic beta cells under endoplasmic reticulum stress.

Thiol-based post-translational modifications (PTMs) play a key role in redox-dependent regulation and signaling. Functional cysteine (Cys) sites serve as redox switches, regulated through multiple types of PTMs. Herein, we aim to characterize the complexity of thiol PTMs at the proteome level through the establishment of a direct detection workflow. The LC-MS/MS based workflow allows for simultaneous quantification of protein abundances and multiple types of thiol PTMs. To demonstrate its utility, the workflow was applied to mouse pancreatic ß-cells (ß-TC-6) treated with thapsigargin to induce endoplasmic reticulum (ER) stress. This resulted in the quantification of >9000 proteins and multiple types of thiol PTMs, including intra-peptide disulfide (S-S), S-glutathionylation (SSG), S-sulfinylation (SO2H), S-sulfonylation (SO3H), S-persulfidation (SSH), and S-trisulfidation (SSSH). Proteins with significant changes in abundance were observed to be involved in canonical pathways such as autophagy, unfolded protein response, protein ubiquitination pathway, and EIF2 signaling. Moreover, ~500 Cys sites were observed with one or multiple types of PTMs with SSH and S-S as the predominant types of modifications. In many cases, significant changes in the levels of different PTMs were observed on various enzymes and their active sites, while their protein abundance exhibited little change. These results provide evidence of independent translational and post-translational regulation of enzyme activity. The observed complexity of thiol modifications on the same Cys residues illustrates the challenge in the characterization and interpretation of protein thiol modifications and their functional regulation.

New Insights Into Lignification via Network and Multi-Omics Analyses of Arogenate Dehydratase Knock-Out Mutants in Arabidopsis thaliana.

Multiple Arabidopsis arogenate dehydratase (ADT) knock-out (KO) mutants, with phenotypes having variable lignin levels (up to circa 70% reduction), were studied to investigate how differential reductions in ADTs perturb its overall plant systems biology. Integrated "omics" analyses (metabolome, transcriptome, and proteome) of wild type (WT), single and multiple ADT KO lines were conducted. Transcriptome and proteome data were collapsed into gene ortholog (GO) data, with this allowing for enzymatic reaction and metabolome cross-comparisons to uncover dominant or likely metabolic biosynthesis reactions affected. Network analysis of enzymes-highly correlated to stem lignin levels-deduced the involvement of novel putative lignin related proteins or processes. These included those associated with ribosomes, the spliceosome, mRNA transport, aminoacyl tRNA biosynthesis, and phosphorylation. While prior work helped explain lignin biosynthesis regulation at the transcriptional level, our data here provide support for a new hypothesis that there are additional post-transcriptional and translational level processes that need to be considered. These findings are anticipated to lead to development of more accurate depictions of lignin/phenylpropanoid biosynthesis models in situ, with new protein targets identified for further biochemical analysis and/or plant bioengineering. Additionally, using KEGG defined functional categorization of proteomics and transcriptomics analyses, we detected significant changes to glucosinolate, α-linolenic acid, nitrogen, carotenoid, aromatic amino acid, phenylpropanoid, and photosynthesis-related metabolic pathways in ADT KO mutants. Metabolomics results also revealed that putative carotenoid and galactolipid levels were generally increased in amount, whereas many glucosinolates and phenylpropanoids (including flavonoids and lignans) were decreased in the KO mutants.

A resource of lipidomics and metabolomics data from individuals with undiagnosed diseases.

Every year individuals experience symptoms that remain undiagnosed by healthcare providers. In the United States, these rare diseases are defined as a condition that affects fewer than 200,000 individuals. However, there are an estimated 7000 rare diseases, and there are an estimated 25-30 million Americans in total (7.6-9.2% of the population as of 2018) affected by such disorders. The NIH Common Fund Undiagnosed Diseases Network (UDN) seeks to provide diagnoses for individuals with undiagnosed disease. Mass spectrometry-based metabolomics and lipidomics analyses could advance the collective understanding of individual symptoms and advance diagnoses for individuals with heretofore undiagnosed disease. Here, we report the mass spectrometry-based metabolomics and lipidomics analyses of blood plasma, urine, and cerebrospinal fluid from 148 patients within the UDN and their families, as well as from a reference population of over 100 individuals with no known metabolic diseases. The raw and processed data are available to the research community so that they might be useful in the diagnoses of current or future patients suffering from undiagnosed disorders.

Dynamic Time-Warping Correction for Shifts in Ultrahigh Resolving Power Ion Mobility Spectrometry and Structures for Lossless Ion Manipulations.

Detection of arrival time shifts between ion mobility spectrometry (IMS) separations can limit achievable resolving power (Rp), particularly when multiple separations are summed or averaged, as commonly practiced in IMS. Such variations can be apparent in higher Rp measurements and are particularly evident in long path length traveling wave structures for lossless ion manipulations (SLIM) IMS due to their typically much longer separation times. Here, we explore data processing approaches employing single value alignment (SVA) and nonlinear dynamic time warping (DTW) to correct for variations between IMS separations, such as due to pressure fluctuations, to enable more effective spectrum summation for improving Rp and detection of low-intensity species. For multipass SLIM IMS separations, where narrow mobility range measurements have arrival times that can extend to several seconds, the SVA approach effectively corrected for such variations and significantly improved Rp for summed separations. However, SVA was much less effective for broad mobility range separations, such as obtained with multilevel SLIM IMS. Changes in ions' arrival times were observed to be correlated with small pressure changes, with approximately 0.6% relative arrival time shifts being common, sufficient to result in a loss of Rp for summed separations. Comparison of the approaches showed that DTW alignment performed similarly to SVA when used over a narrow mobility range but was significantly better (providing narrower peaks and higher signal intensities) for wide mobility range data. We found that the DTW approach increased Rp by as much as 115% for measurements in which 50 IMS separations over 2 s were summed. We conclude that DTW is superior to SVA for ultra-high-resolution broad mobility range SLIM IMS separations and leads to a large improvement in effective Rp, correcting for ion arrival time shifts regardless of the cause, as well as improving the detectability of low-abundance species. Our tool is publicly available for use with universal ion mobility format (.UIMF) and text (.txt) files.

Proteogenomic and metabolomic characterization of human glioblastoma.

Glioblastoma (GBM) is the most aggressive nervous system cancer. Understanding its molecular pathogenesis is crucial to improving diagnosis and treatment. Integrated analysis of genomic, proteomic, post-translational modification and metabolomic data on 99 treatment-naive GBMs provides insights to GBM biology. We identify key phosphorylation events (e.g., phosphorylated PTPN11 and PLCG1) as potential switches mediating oncogenic pathway activation, as well as potential targets for EGFR-, TP53-, and RB1-altered tumors. Immune subtypes with distinct immune cell types are discovered using bulk omics methodologies, validated by snRNA-seq, and correlated with specific expression and histone acetylation patterns. Histone H2B acetylation in classical-like and immune-low GBM is driven largely by BRDs, CREBBP, and EP300. Integrated metabolomic and proteomic data identify specific lipid distributions across subtypes and distinct global metabolic changes in IDH-mutated tumors. This work highlights biological relationships that could contribute to stratification of GBM patients for more effective treatment.

Unveiling molecular signatures of preeclampsia and gestational diabetes mellitus with multi-omics and innovative cheminformatics visualization tools.

To fully enable the development of diagnostic tools and progressive pharmaceutical drugs, it is imperative to understand the molecular changes occurring before and during disease onset and progression. Systems biology assessments utilizing multi-omic analyses (e.g. the combination of proteomics, lipidomics, genomics, etc.) have shown enormous value in determining molecules prevalent in diseases and their associated mechanisms. Herein, we utilized multi-omic evaluations, multi-dimensional analysis methods, and new cheminformatics-based visualization tools to provide an in depth understanding of the molecular changes taking place in preeclampsia (PRE) and gestational diabetes mellitus (GDM) patients. Since PRE and GDM are two prevalent pregnancy complications that result in adverse health effects for both the mother and fetus during pregnancy and later in life, a better understanding of each is essential. The multi-omic evaluations performed here provide new insight into the end-stage molecular profiles of each disease, thereby supplying information potentially crucial for earlier diagnosis and treatments.

Correction: Proteogenomic Characterization of Ovarian HGSC Implicates Mitotic Kinases, Replication Stress in Observed Chromosomal Instability.

[This corrects the article PMC7289043.].

Proteogenomic Characterization of Ovarian HGSC Implicates Mitotic Kinases, Replication Stress in Observed Chromosomal Instability.

In the absence of a dominant driving mutation other than uniformly present TP53 mutations, deeper understanding of the biology driving ovarian high-grade serous cancer (HGSC) requires analysis at a functional level, including post-translational modifications. Comprehensive proteogenomic and phosphoproteomic characterization of 83 prospectively collected ovarian HGSC and appropriate normal precursor tissue samples (fallopian tube) under strict control of ischemia time reveals pathways that significantly differentiate between HGSC and relevant normal tissues in the context of homologous repair deficiency (HRD) status. In addition to confirming key features of HGSC from previous studies, including a potential survival-associated signature and histone acetylation as a marker of HRD, deep phosphoproteomics provides insights regarding the potential role of proliferation-induced replication stress in promoting the characteristic chromosomal instability of HGSC and suggests potential therapeutic targets for use in precision medicine trials.

Block Design with Common Reference Samples Enables Robust Large-Scale Label-Free Quantitative Proteome Profiling.

Label-free quantitative proteomics has become an increasingly popular tool for profiling global protein abundances. However, one major limitation is the potential performance drift of the LC-MS platform over time, which, in turn, limits its utility for analyzing large-scale sample sets. To address this, we introduce an experimental and data analysis scheme based on a block design with common references within each block for enabling large-scale label-free quantification. In this scheme, a large number of samples (e.g., >100 samples) are analyzed in smaller and more manageable blocks, minimizing instrument drift and variability within individual blocks. Each designated block also contains common reference samples (e.g., controls) for normalization across all blocks. We demonstrated the robustness of this approach by profiling the proteome response of human macrophage THP-1 cells to 11 engineered nanomaterials at two different doses. A total of 116 samples were analyzed in six blocks, yielding an average coverage of 4500 proteins per sample. Following a common reference-based correction, 2537 proteins were quantified with high reproducibility without any imputation of missing values from 116 data sets. The data revealed the consistent quantification of proteins across all six blocks, as illustrated by the highly consistent abundances of house-keeping proteins in all samples and the high levels of correlation among samples from different blocks. The data also demonstrated that label-free quantification is robust and accurate enough to quantify even very subtle abundance changes as well as large fold-changes. Our streamlined workflow is easy to implement and can be readily adapted to other large cohort studies for reproducible label-free proteome quantification.

Accurate Identification of Deamidation and Citrullination from Global Shotgun Proteomics Data Using a Dual-Search Delta Score Strategy.

Proteins with deamidated/citrullinated amino acids play critical roles in the pathogenesis of many human diseases; however, identifying these modifications in complex biological samples has been an ongoing challenge. Herein we present a method to accurately identify these modifications from shotgun proteomics data generated by a deep proteome profiling study of human pancreatic islets obtained by laser capture microdissection. All MS/MS spectra were searched twice using MSGF+ database matching, with and without a dynamic +0.9840 Da mass shift modification on amino acids asparagine, glutamine, and arginine (NQR). Consequently, each spectrum generates two peptide-to-spectrum matches (PSMs) with MSGF+ scores, which were used for the Delta Score calculation. It was observed that all PSMs with positive Delta Score values were clustered with mass errors around 0 ppm, while PSMs with negative Delta Score values were distributed nearly equally within the defined mass error range (20 ppm) for database searching. To estimate false discovery rate (FDR) of modified peptides, a "target-mock" strategy was applied in which data sets were searched against a concatenated database containing "real-modified" (+0.9840 Da) and "mock-modified" (+1.0227 Da) peptide masses. The FDR was controlled to ∼2% using a Delta Score filter value greater than zero. Manual inspection of spectra showed that PSMs with positive Delta Score values contained deamidated/citrullinated fragments in their MS/MS spectra. Many citrullinated sites identified in this study were biochemically confirmed as autoimmunogenic epitopes of autoimmune diseases in literature. The results demonstrated that in situ deamidated/citrullinated peptides can be accurately identified from shotgun tissue proteomics data using this dual-search Delta Score strategy. Raw MS data is available at ProteomeXchange (PXD010150).

Proteogenomic Characterization of Endometrial Carcinoma.

We undertook a comprehensive proteogenomic characterization of 95 prospectively collected endometrial carcinomas, comprising 83 endometrioid and 12 serous tumors. This analysis revealed possible new consequences of perturbations to the p53 and Wnt/ß-catenin pathways, identified a potential role for circRNAs in the epithelial-mesenchymal transition, and provided new information about proteomic markers of clinical and genomic tumor subgroups, including relationships to known druggable pathways. An extensive genome-wide acetylation survey yielded insights into regulatory mechanisms linking Wnt signaling and histone acetylation. We also characterized aspects of the tumor immune landscape, including immunogenic alterations, neoantigens, common cancer/testis antigens, and the immune microenvironment, all of which can inform immunotherapy decisions. Collectively, our multi-omic analyses provide a valuable resource for researchers and clinicians, identify new molecular associations of potential mechanistic significance in the development of endometrial cancers, and suggest novel approaches for identifying potential therapeutic targets.

From Prevention to Disease Perturbations: A Multi-Omic Assessment of Exercise and Myocardial Infarctions.

While a molecular assessment of the perturbations and injury arising from diseases is essential in their diagnosis and treatment, understanding changes due to preventative strategies is also imperative. Currently, complex diseases such as cardiovascular disease (CVD), the leading cause of death worldwide, suffer from a limited understanding of how the molecular mechanisms taking place following preventive measures (e.g., exercise) differ from changes occurring due to the injuries caused from the disease (e.g., myocardial infarction (MI)). Therefore, this manuscript assesses lipidomic changes before and one hour after exercise treadmill testing (ETT) and before and one hour after a planned myocardial infarction (PMI) in two separate patient cohorts. Strikingly, unique lipidomic perturbations were observed for these events, as could be expected from their vastly different stresses on the body. The lipidomic results were then combined with previously published metabolomic characterizations of the same patients. This integration provides complementary insights into the exercise and PMI events, thereby giving a more holistic understanding of the molecular changes associated with each.

SLIM Ultrahigh Resolution Ion Mobility Spectrometry Separations of Isotopologues and Isotopomers Reveal Mobility Shifts due to Mass Distribution Changes.

We report on separations of ion isotopologues and isotopomers using ultrahigh-resolution traveling wave-based Structures for Lossless Ion Manipulations with serpentine ultralong path and extended routing ion mobility spectrometry coupled to mass spectrometry (SLIM SUPER IMS-MS). Mobility separations of ions from the naturally occurring ion isotopic envelopes (e.g., [M], [M+1], [M+2], ... ions) showed the first and second isotopic peaks (i.e., [M+1] and [M+2]) for various tetraalkylammonium ions could be resolved from their respective monoisotopic ion peak ([M]) after SLIM SUPER IMS with resolving powers of ∼400-600. Similar separations were obtained for other compounds (e.g., tetrapeptide ions). Greater separation was obtained using argon versus helium drift gas, as expected from the greater reduced mass contribution to ion mobility described by the Mason-Schamp relationship. To more directly explore the role of isotopic substitutions, we studied a mixture of specific isotopically substituted (15N, 13C, and 2H) protonated arginine isotopologues. While the separations in nitrogen were primarily due to their reduced mass differences, similar to the naturally occurring isotopologues, their separations in helium, where higher resolving powers could also be achieved, revealed distinct additional relative mobility shifts. These shifts appeared correlated, after correction for the reduced mass contribution, with changes in the ion center of mass due to the different locations of heavy atom substitutions. The origin of these apparent mass distribution-induced mobility shifts was then further explored using a mixture of Iodoacetyl Tandem Mass Tag (iodoTMT) isotopomers (i.e., each having the same exact mass, but with different isotopic substitution sites). Again, the observed mobility shifts appeared correlated with changes in the ion center of mass leading to multiple monoisotopic mobilities being observed for some isotopomers (up to a ∼0.04% difference in mobility). These mobility shifts thus appear to reflect details of the ion structure, derived from the changes due to ion rotation impacting collision frequency or momentum transfer, and highlight the potential for new approaches for ion structural characterization.

Lipid Mini-On: mining and ontology tool for enrichment analysis of lipidomic data.

SUMMARY: Here we introduce Lipid Mini-On, an open-source tool that performs lipid enrichment analyses and visualizations of lipidomics data. Lipid Mini-On uses a text-mining process to bin individual lipid names into multiple lipid ontology groups based on the classification (e.g. LipidMaps) and other characteristics, such as chain length. Lipid Mini-On provides users with the capability to conduct enrichment analysis of the lipid ontology terms using a Shiny app with options of five statistical approaches. Lipid classes can be added to customize the user's database and remain updated as new lipid classes are discovered. Visualization of results is available for all classification options (e.g. lipid subclass and individual fatty acid chains). Results are also visualized through an editable network of relationships between the individual lipids and their associated lipid ontology terms. The utility of the tool is demonstrated using biological (e.g. human lung endothelial cells) and environmental (e.g. peat soil) samples. AVAILABILITY AND IMPLEMENTATION: Rodin (R package: https://github.com/PNNL-Comp-Mass-Spec/Rodin), Lipid Mini-On Shiny app (https://github.com/PNNL-Comp-Mass-Spec/LipidMiniOn) and Lipid Mini-On online tool (https://omicstools.pnnl.gov/shiny/lipid-mini-on/). SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.